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ABSTRACT 


The  signal  classification  performance  of  four  types  of  Electronics  Support  Mea¬ 
sure  (ESM)  communications  detection  systems  is  compared  from  the  standpoint  of 
the  unintended  receiver  (interceptor).  Typical  digital  communication  signals  con¬ 
sidered  include  Binary  Phase  Shift  Keying  (BPSK),  Quadrature  Phase  Shift  Keying 
(QPSK),  Frequency  Shift  Keying  (FSK),  and  On-Off  Keying  (OOK).  The  analysis 
emphasizes  the  use  of  available  signal  processing  software.  Detection  methods  com¬ 
pared  include  broadband  energy  detection,  FFT-based  narrowband  energy  detection, 
and  two  correlation  methods  which  employ  the  Fast  Fourier  Transform  (FFT).  The 
correlation  methods  utilize  modified  time-frequency  distributions,  where  one  of  these 
is  based  on  the  Wigner-Ville  Distribution  (WVD).  Gaussian  white  noise  is  added  to 
the  signal  to  simulate  various  signal-to- noise  ratios  (SNRs). 
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I.  INTRODUCTION 


Communication  signal  detection  and  interception  is  an  important  aspect  of  mod¬ 
ern  Command  and  Control  Warfare  (C2W).  Detection  and  interception  techniques 
extract  signal  features  such  as  the  carrier  frequency,  bit  rate,  type  of  modulation,  and 
chip  rate  in  the  case  of  spread  spectrum  signals  [1,  2]. 

The  detection  problem  for  the  unintended  interceptor  or  receiver  is  more  difficult 
than  that  for  the  intended  receiver,  since  the  interceptor  may  not  know  the  carrier 
frequency,  the  type  of  modulation,  or  the  message  bit  transition  instants.  The  work 
done  in  this  thesis  provides  some  basic  idea^  on  extracting  unknown  signal  parameters 
or  features  such  as  modulation  type,  bandwidth,  bit  rate,  and  carrier  frequency. 

Four  types  of  detectors  are  implemented:  the  broadband  energy  detector  with 
adjustable  integration  time,  the  FFT-based  detector  with  adjustable  integration  time, 
and  two  types  of  correlation  detectors  based  on  a  modified  time-frequency  distribution 
method.  The  latter  two  are  denoted  as  method-1,  which  is  related  to  the  Fourier 
transform  of  lagged  products  (over  lag  m),  and  method-2  which  is  based  on  the 
Fourier  transform  of  lagged  products  (over  time  n). 

These  four  signal  detectors  are  compared  with  respect  to  their  ability  to  dis¬ 
tinguish  among  four  types  of  modulations  (BPSK,  QPSK,  FSK,  OOK).  The  overall 
schematic  for  the  systems  is  shown  in  Figure  1. 

Chapter  II  presents  the  signal  generation  and  the  underlying  assumptions. 
Chapter  III  describes  the  detectors.  Chapter  IV  presents  experiments  of  the  de¬ 
tection  methods  for  each  modulation  type  as  a  function  of  SNR.  The  performance 
of  the  four  detectors  is  summarized  and  recommendations  for  future  studies  are  dis¬ 
cussed  in  Chapter  V.  The  MATLAB  programs  used  for  this  study  are  listed  in  the 
Appendix. 
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Interception  of  spread  spectrum  signals  is  not  included  in  this  thesis  which  is 
limited  to  detectors  for  signals  with  positive  SNR.  The  intercept  problem  is  partic¬ 
ularly  difficult  for  spread  spectrum  signals  which  typically  have  negative  SNRs.  In 
this  case  techniques  such  as  “prefilter  delay  and  multiply”  introduce  nonlinear  trans¬ 
formations  which  enable  certain  feature  extractions  [3] . 
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II.  SIGNAL  GENERATION 


A.  INTRODUCTION 

In  digital  communications  systems,  various  types  of  modulation  techniques  are 
employed  to  transmit  information.  Some  of  the  more  common  forms  of  modulation 
are  amplitude  modulation  (AM),  frequency  modulation  (FM),  and  phase  modulation 
(PM).  In  AM,  FM,  or  PM,  the  amplitude,  frequency,  or  phase,  respectively,  is  varied 
by  the  modulating  signal  [4]. 

B.  SIGNAL  GENERATION 

Four  types  of  binary  modulation  schemes  (BPSK,  OOK,  FSK,  QPSK)  are  used 
in  the  work  addressed  by  this  thesis.  All  messages  have  the  same  sine  wave  carrier 
and  a  peak  amplitude  of  one.  A  common  16-bit  information  sequence  is  input  to 
each  modulator,  resulting  in  four  carrier  cycles  per  bit.  There  are  four  sampling 
points  per  cycle,  resulting  in  a  “message”  length  of  256  points.  These  parameters 
remain  constant  throughout  the  simulations.  To  allow  simulation  of  corrupted  data, 
Gaussian  noise  is  added  to  create  10-decibel  (dB)  and  0-dB  SNRs. 

1.  Binary  Phase  Shift  Keying  (BPSK) 

This  signal  is  a  common  example  of  phase  modulation,  where  the  phase  of  the 
carrier  is  shifted  180  degrees  relative  to  the  quiescent  carrier  phase  during  the  zero  of 
the  modulating  signal.  In  Figure  2  the  binary  values  of  -1-1  and  0  are  represented  by 
-|-1  and  —1,  respectively.  Figure  2  shows  the  modulating  signal  and  resulting  BPSK 
signal  while  Figure  3a  shows  the  BPSK  signal  at  an  SNR  of  10-dB.  Note  that  it  is 
diflScult  to  estimate  the  type  of  modulation.  Figure  3b  shows  the  BPSK  signal  at  an 
SNR  of  0-dB.  Note  that  it  is  impossible  to  separate  the  signal  from  the  noise  . 
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2.  On-OfF  Keying  (OOK) 

This  type  of  modulation  technique  is  also  called  amplitude  shift  keying  (ASK). 
It  consists  of  keying  (switching)  a  carrier  sinusoid  ON  and  OFF  with  a  unipolar  binary 
signal.  It  is  identical  to  unipolar  binary  amplitude  modulation  resulting  in  a  Double 
Sideband-Suppressed  Carrier  (DSB-SC)  signal.  Figure  4  shows  the  modulating  signal 
and  the  resulting  OOK  signal,  while  Figure  5a  shows  the  OOK  signal  at  an  SNR  of 
10-dB.  In  this  figure  the  features  of  OOK  modulation  can  be  recognized.  In  contrast, 
as  shown  in  Figure  5b,  at  0-dB  SNR,  it  is  difficult  to  differentiate  between  the  signal 
and  noise. 

3.  Quadrature  Phase  Shift  Keying  (QPSK) 

This  signal  can  be  considered  as  two  BPSK  signals  occupying  the  same  channel, 
in  a  quadrature  arrangement.  The  QPSK  signals  are  created  by  90  degrees  relative 
phase  shifts,  as  shown  in  Figure  6.  The  QPSK  signal  at  an  SNR  of  10-dB  is  shown 
in  Figure  7a.  The  signal  can  be  recognized.  Figure  7b  shows  the  QPSK  signal  at  an 
SNR  of  0-dB,  which  does  not  allow  differentiation  of  the  signal  from  the  noise. 

4.  Frequency  Shift  Keying  (FSK) 

In  this  case,  the  frequency  of  a  carrier  will  shift  from  a  mark  frequency  (corre¬ 
sponding  to  a  binary  1)  to  a  space  frequency  (corresponding  to  a  binary  0)  according 
to  the  baseband  modulation  signal.  This  is  illustrated  in  Figure  8,  which  shows  the 
modulation  signal  and  the  resulting  FSK  signal.  In  the  experiments,  the  difference 
between  the  mark  frequency  and  space  frequency  is  one  carrier  cycle  per  bit  (i.e.,  the 
mark  frequency  is  equal  to  5  cycles  and  the  space  frequency  is  equal  to  4  cycles  per 
bit).  The  FSK  signal  at  an  SNR  of  10-dB,  shown  in  Figure  9a,  can  still  be  recognized. 
Figure  9b  shows  the  FSK  at  an  SNR  of  0-dB,  which  does  not  allow  differentiation  of 
the  signal  from  the  noise. 
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III.  SIGNAL  INTERCEPT  DETECTION 

A.  INTRODUCTION 

In  the  detection  of  signals,  especially  LPI  (low  probability  of  intercept)  signals, 
the  ease  of  detection  and  the  methods  used  vary  greatly,  depending  upon  whether 
the  receiver  is  the  intended  or  unintended  receiver.  For  the  purpose  of  this  thesis,  the 
intended  receiver  is  the  person  or  system  that  receives  the  signal  from  a  cooperative 
transmitter  and  has  full  access  to  the  frequency  and  code  information  required  for 
detection  and  recovery  of  the  message.  The  unintended  receiver,  on  the  other  hand, 
does  not  have  access  to  the  frequency  and  code  information  and  is  attempting  to 
detect  the  presence  of  the  signal  and/or  recover  the  message  contained  in  the  data 
[5]. 

The  traditional  method  of  signal  intercept  detection  employed  by  an  unintended 
receiver  is  the  radiometer  (integration  of  energy)  detector  which  requires  a  positive 
SNR.  In  recent  years  this  has  been  extended  to  intercept  spread  spectrum  signals 
with  negative  SNR  [3].  Negative  SNR  requires  a  long  observation  time  and  accurate 
estimation  of  the  threshold  set  by  the  noise  above.  Only  positive  SNR  detection  is 
considered  in  this  work.  Four  types  of  detectors  are  described  in  the  following  section: 
energy,  FFT,  method-1,  and  method-2.  The  broadband  energy  detector  is  intended 
to  be  a  benchmark  for  comparison  with  the  other  detectors  and  is  used  for  OOK  only. 
The  FFT  detector  is  more  effective  than  the  energy  detector.  Methods  1  and  2  are 
based  on  time-frequency  distribution.  One  of  these  methods  is  related  to  the  WVD. 
The  various  detectors  are  illustrated  in  Figure  1. 
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B.  BROADBAND  ENERGY  DETECTOR 


The  interceptor  employs  the  radiometer  or  traditional  energy  detector  to  deter¬ 
mine  the  presence  of  a  signal  which  may  have  an  unknown  bit  rate  or  start  time.  The 

energy  detector  integrates  the  signal  energy  for  a  finite  length  of  time  as  described  by 

I+N-l 

E=  T.  WOP.  (1) 

«=/ 

where  E  is  the  energy  and  r(i)  is  the  received  signal  at  time  i.  The  summation  is 
over  N  points,  where  N  is  the  segment  length,  and  I  is  the  starting  point.  If  the  total 
energy,  received  over  a  given  integration  length  N,  is  greater  than  some  threshold 
value  which  is  a  function  of  ambient  noise,  then  a  “detection”  occurs.  The  threshold 
depends  on  the  acceptable  false  alarm  rate  and  on  the  probability  of  detection.  When 
the  noise  is  non-stationary  or  strong  interference  signals  (ISs)  are  present,  frequently 
the  case  in  practice,  this  method  is  ineffective  [4,  7,  8]. 

The  energy  detector  described  above  is  used  on  the  OOK  signal  modulation 
type.  Some  detection  examples  are  shown  in  Chapter  IV  as  a  function  of  integration 
length  A,  start  time  /,  and  SNR. 

C.  FAST  FOURIER  TRANSFORM  (FFT)  DETECTION 

The  Fast  Fourier  Transform  (FFT)  algorithm  is  used  in  the  analysis  of  finite 
duration  sequences.  In  this  thesis,  an  attempt  is  made  to  use  the  FFT  to  extract 
parameters  such  as  modulation  type,  bit  rate,  transition  time,  band  width,  etc.,  from 
the  intercepted  signals  [7].  This  method  is  more  efficient  than  broadband  energy  de¬ 
tection.  The  FFT  parameters  are  selected  by  the  user  to  enhance  the  signal  frequency 
components  by  using  a  range  of  likely  integration  times.  This  method  is  equivalent  to 
a  parallel  bank  of  noncoherent  correlation  detectors  in  which  the  detection  decision 
is  based  on  comparing  the  detector  outputs  to  the  threshold.  This  is  similar  to  the 
method  used  to  noncoherently  detect  and  demodulate  M-ary  Frequency  Shift  Keying 
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(MFSK)  [2,  4,  8,  9].  The  algorithm  used  for  calculating  the  detector  energy  output 

level  E{k)  at  each  frequency  bin  k  is 

Nt  I+N-1 

E{k)  =  £  E 

i=i 

where  i  is  the  time  index,  E{k)  is  the  total  output  energy  in  the  time  interval  and 
frequency  interval  with  index  k,  r(i)  is  the  received  signal  at  time  i,  I  corresponds 
to  the  integration  or  observation  start  time,  Nt  is  the  number  of  terms  used  in  the 
power  average,  and  N  is  the  coherent  integration  length.  The  FFT  intervals  (i.e., 
(/,  /  +  A^))  are  either  contiguous  or  overlapping. 

D.  MODIFIED  TIME-FREQUENCY  DISTRIBUTION  METHODS 

One  remarkable  fact  regarding  time-frequency  distribution  is  that  many  deriva¬ 
tions  and  approaches  have  been  suggested.  It  is  important  to  understand  the  ideas  and 
arguments  that  have  been  given,  as  variations  and  insights  into  them  will  undoubt¬ 
edly  lead  the  way  to  further  development.  These  approaches  will  not  be  presented  in 
this  work,  but  we  will  refer  to  them  for  further  study  [10]. 

Two  types  of  time-frequency  distributions  are  presented  in  this  work,  method- 
1,  which  is  related  to  the  Fourier  transform  of  lagged  products  (over  lag  m),  and 
method-2,  which  is  based  on  the  Fourier  transform  of  lagged  products  (over  time 
n).  The  output  magnitude  of  each  method  is  examined  to  evaluate  its  diagnostic 
potential  in  recognizing  the  modulation  of  the  signal. 

1.  Method-1 
Method- 1  uses 

N-x 

X{k,n)  =  ^  x{n)x{n  -f  ,  (3) 

m=0 

where  n  is  independent  of  delay  time  m.  Eq.  (3)  is  interpreted  as 

N-l 

X{k,  n)  =  x{n)  ^  x(n  +  ,  (4) 

m=0 
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where  k  is  the  frequency  index,  n  the  time  index,  m  the  time  delay,  and  N  the  number 
of  samples  during  the  observation  period.  x(n)  is  the  estimated  value  of  the  mean  of 
x{n)  simply  given  by 

_  1  n-\‘N—\ 

=  Y.  ^(*)-  (5) 

i=:n 

2.  Method-2 

The  following  algorithm  is  used  to  differentiate  between  the  modulated  signals 

L— m— 1 

X{k,m)  =  ^  x(n)a;(n  +  ,  (6) 

n=0 

where  k  is  the  frequency  index,  n  the  time  index,  m  is  the  delay,  and  L  the  number 
of  samples  used  during  the  observation  period.  Even  though  a  spectrum  can  be  com¬ 
puted  for  every  m  (m  =  0, 1,  •  •  •),  we  choose  to  compute  for  computational  efficiency, 
the  spectra  only  at  evenly  spaced  points.  That  is,  m  might  be  indexed  to  go  from  mo, 
to  mo  -hstep,  to  mo  -t-2  step,  and  so  on.  The  resulting  time-frequency  surface  will  be 
computed  much  faster,  yet  not  too  much  information  is  lost  if  the  nonstationarity  of 
the  signal  and  the  step  size  are  compatible. 

As  mentioned  earlier,  this  method  is  related  to  WVD,  which  is  a  three-dimensional 
(energy  vs.  time  and  frequency)  representation  of  a  signal  that  is  particularly  well- 
suited  for  analysis  of  non-stationary  signals.  It  was  first  introduced  by  Eugene  Wigner 
and  then  applied  to  signal  analysis  in  1948  by  J.  Ville. 

The  WVD  is  a  transformation  of  a  continuous  time  signal  into  the  time-frequency 
domain.  The  auto  WVD  is  defined  as  [6,  10] 

WVDrAt,'^)  =  /_”  ‘-’“'r  ((  +  j)  r-  (i  -^)dT.  (7) 

In  terms  of  the  transform  variables,  it  is 

(“  +  0  (8) 
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where 


r{t)  =  complex  valued  signal  in  time  domain 
R{u!)  =  complex  valued  signal  in  frequency  domain 

t  =  time 

(jj  =  frequency 

*  =  complex  conjugate 

Since  the  concern  in  this  work  is  to  detect  and  extract  parameters  from  signals, 
there  are  properties  of  WVD  which  should  be  kept  in  mind  when  interpreting  the 
processing  results  [10]: 

(1)  WVD  always  goes  to  zero  at  the  beginning  and  the  end  of  a  finite-duration 
signal,  sometimes  called  the  support  properties  of  WVD. 

(2)  In  general,  the  Wigner  distribution  is  not  zero  when  the  signal  is  zero. 
This  causes  considerable  difficulty  in  interpretation. 

(3)  According  to  Cohen  [10],  “In  general,  the  Wigner  distribution  has  a  noisy 
behavior,  that  is,  if  there  is  a  noise  in  a  small  finite  time  of  the  signal,  this 
noise  will  appear  for  all  time.” 
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IV.  EXPERIMENTAL  PERFORMANCE  RESULTS 


A.  BASIC  CONCEPTS 

The  energy  detector  is  only  tested  using  the  OOK  signal  modulation.  The  other 
three  detection  methods  are  tested  by  comparing  the  detector’s  three-dimensional 
output  display.  The  ability  of  a  human  to  determine  signal  presence,  modulation 
type,  and  other  signal  features  from  the  displays  provides  useful  insights  for  future 
design  of  automatic  digital  classification  algorithms.  The  experiment  is  conducted 
as  illustrated  in  Figure  1.  At  the  left-hand  side,  the  signals  (BPSK,  OOK,  QPSK, 
FSK)  are  generated  and  saved  in  files  for  input  to  the  various  detectors  that  are  to 
be  compared.  This  assures  a  fixed  input  waveform  as  the  detector  parameters  (e.g., 
FFT  window  size)  are  varied. 

The  broadband  energy  detector  is  implemented  with  the  ENERGY. m  program  and 
the  FFT  detector  with  SHIFT2.m.  The  time-frequency  distribution  method- 1  and 
method-2  are  implemented  with  the  CORRl  .m  and  C0RR2.m  programs,  respectively.  In 
all  plots  the  frequency  axis  is  calibrated  in  FFT  bin  numbers. 

B.  BROADBAND  ENERGY  DETECTION 

As  mentioned  earlier,  we  restrict  our  broadband  energy  detection  experiments  to 
the  OOK  signal  because  the  other  modulation  techniques  have  continuous  time  energy, 
while  the  OOK  signal  has  an  ON-OFF  energy  representation.  This  mahes  OOK 
message  identification  simple,  especially  at  high  SNR.  Due  to  the  time  constraints, 
we  did  not  compare  the  outputs  with  signal  ON  and  signal  OFF. 

For  analysis  purposes,  the  output  of  the  energy  detector,  with  different  noise 
levels,  will  be  represented  by  two  plot  types.  First,  an  overall  plot,  shown  in  Figure 
10a,  11a,  and  12a  displays  the  detector  output  as  a  function  of  integration  time 
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(samples)  and  time.  The  second  plot  in  these  figures  shows  the  detector  output  for 
selected  integration  time  lengths  superimposed  on  each  other.  This  will  demonstrate 
the  effect  of  the  integration  length  on  the  detector  output,  as  shown  in  Figures  10b, 
lib,  and  12b.  For  comparison  purposes,  the  detector  input  signal  is  shown  in  Figures 
10c,  11c,  and  12c. 

Figure  10a  shows  the  energy  detector  output  for  the  noise- free  OOK  signal. 
Note  that  the  shorter  the  integration  time,  the  steeper  become  the  leading  and  the 
traihng  edges  of  the  output,  approaching  a  rectangular  shape.  Distinct  bits  result  in 
a  triangular  output  shape  when  the  integration  time  length  is  exactly  equal  to  the 
bit  length.  To  clarify  this  idea,  four  integration  lengths  are  selected  from  Figure  10a 
and  superimposed  in  Figure  10b. 

Comparing  Figures  10a  and  10b,  it  is  noted  that  the  amplitude  increases  and 
the  output  starts  earlier  with  increasing  integration  time.  The  amplitude  increa.se  is 
caused  by  summing  more  non-negative  numbers  as  the  integration  time  is  increased. 

Comparing  the  input  and  the  output  of  the  energy  detector,  note  that,  for  a 
short  integration  length,  the  output  of  the  energy  detector  matches  the  amplitude 
modulated  signal.  Hence,  under  noise-free  conditions,  the  pulse  width  and  the  bit 
rate  can  readily  be  estimated. 

The  output  of  the  energy  detector  for  an  SNR  of  10-dB  as  a  function  of  different 
integration  times  is  shown  in  Figure  11a.  Note  that  it  is  difficult  to  detect  or  separate 
the  signal  from  the  noise  at  short  integration  lengths  (integration  lengths  1  through 
4).  Using  more  samples  (i.e.,  longer  integration  time,  say  6  through  16)  makes  the 
signal  recognizable. 

Some  results  for  different  integration  times  are  shown  in  Figure  11b.  Note  that 
at  lengths  8  and  16,  the  modulation  signal  can  be  recognized.  In  contrast,  when  using 
2  or  4  samples,  the  signal  cannot  be  detected. 
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As  shown  in  Figures  12a  and  12b,  the  energy  detector  output  will  have  a  very 
high  level  of  energy,  which  causes  the  signal  to  be  completely  lost,  as  the  level  of  the 
noise  is  further  increased  (i.e.,  0-dB). 

Figure  12c  shows  the  input  signal  for  comparison  purposes.  It  is  clear  from  this 
that  the  broadband  energy  detector  can  be  used  as  a  pulse-width  detector  at  SNR 
levels  greater  or  equal  to  10-dB. 

C.  FFT  BASED  DETECTION 

BPSK,  OOK,  QPSK,  and  FSK  signals  are  applied  to  the  FFT  detector  to  de¬ 
tect  and  identify  message  features  such  as  the  modulation  type,  carrier  frequency, 
frequency  shift  (if  appropriate),  and  bit  rate. 

All  FFT  detector  simulations  use  an  integration  length  of  16.  This  length  is 
chosen  to  determine  which  message  parameters  can  be  extracted  if  the  FFT  size 
equals  the  number  of  samples  per  message  bit. 

The  user  ha.s  the  option  to  select  an  arbitrary  overlap  factor.  A  75%  overlap 
factor  is  chosen,  hence  a  new  segment  will  start  every  four  points.  The  purpose  of 
using  an  overlap  is  to  enhance  the  signal  frequency  components  along  the  time  axis  to 
extract  more  parameter  information  from  the  signal,  which  might  be  lost  when  using 
a  window. 

Adding  noise  to  the  message  gives  us  a  basic  idea  as  to  whether  the  signal  can 
be  identified  in  the  presence  of  noise.  The  SNR  selection  is  limited  to  two  values, 
10-dB  and  0-dB.  All  messages  are  of  the  same  length,  and  have  the  same  number  of 
points  per  cycle,  number  of  cycles  per  bit,  and  number  of  bits  per  message  length. 

The  output  of  the  FFT  detector  for  the  modulated  signals  is  presented  in  three 
types  of  plots.  The  three-dimensional  plot  of  the  ‘a’  part  of  Figures  13  through  15, 
17  through  19,  21  through  23,  and  25  through  27  represents  the  spectrogram  of  the 
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Amplitude 


signals.  The  ‘b’  parts  of  these  figures  are  the  corresponding  contour  plots.  Figures  16, 
20,  24,  and  28  consist  of  six  subplots.  The  first  column  shows  the  output  averaged  over 
time  and  the  second  column  shows  the  output  averaged  over  frequency  for  different 
levels  of  SNR.  The  time  axis  (samples  or  delay)  is  scaled  by  a  factor  of  4  relative 
to  the  original  time  series.  It  represents  the  shift  of  four  points  between  the  signal 
segments  (i.e.,  75%  overlap  at  a  transform  length  of  16). 

1.  BPSK 

The  output  of  the  detector  driven  by  a  noise-free  BPSK  signal  is  shown  in 
Figure  13a.  A  modulated  carrier  is  indicated.  Spreading  of  the  spectrum  (over 
frequency)  could  indicate  the  presence  of  phase  shift  keying.  Looking  at  Figure  13b, 
which  shows  the  corresponding  contour  plot,  we  note  that  transition  points  (i.e.,  pha^e 
shifts)  are  fairly  easy  to  identify.  A  modulated  carrier  is  present  at  spectral  locations  5 
and  11.  The  bit  rate  can  be  estimated  by  measuring  the  minimum  distance  between 
the  leading  edges  of  two  consecutive  transition  points.  This  minimum  distance  is 
properly  defined  if  all  distances  between  transition  points  (or  the  majority  of  them) 
are  integer  multiples  of  this  minimum  distance.  The  message,  with  a  plus/minus  sign 
uncertainty,  can  be  recovered. 

At  an  SNR  level  of  10-dB,  the  output  of  the  detector  still  shows  the  modulated 
carrier,  as  shown  in  Figure  14a.  Figure  14b  is  a  contour  plot  of  Figure  14a.  It  reveals 
that  the  modulated  carrier  is  at  spectral  locations  5  and  11,  the  transition  points  can 
be  obtained,  and  hence  the  bit  rate  can  still  be  estimated. 

At  an  SNR  level  of  0-dB  the  features  of  the  output  completely  change.  Figure 
15a  shows  the  output  of  the  detector  and  indicates  the  difficulty  in  identifying  the 
output  (i.e.,  any  type  of  modulation  can  cause  it).  Figure  15b  shows  the  corresponding 
contour  plot  of  the  output.  It  appears  that  the  effect  of  noise  at  this  level  may  lead 
us  to  conclude  that  a  FSK  or  OOK  modulated  signal  is  present. 
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Figure  13:  Output  of  FFT  detector  for  BPSK  signal,  no  noise;  75%  overlap 
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Figure  14:  Output  of  FFT  detector  for  BPSK  signal  with  10-dB  SNR;  75%  overlap 
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Figure  16  shows  the  output  averaged  over  frequency  and  time.  At  all  SNR 
levels,  the  carrier  frequency  and  bandwidth  can  be  estimated  using  the  time  averaged 
outputs  (la,  2a,  3a). 

The  frequency  averaged  output  (see  lb,  2b,  and  3b)  provides  little  information. 
Plot  lb  (no  noise)  shows  constant  energy.  At  10-dB  SNR  (see  2b),  it  is  difficult  to 
estimate  the  presence  of  the  BPSK  and  the  level  of  the  noise  completely  masks  the 
signal  at  0-dB  SNR  (see  3b). 

2.  OOK 

This  modulation  technique  is  much  easier  to  identify  than  the  other  types  of 
modulation.  The  output  of  the  FFT  detector  driven  by  a  noise-free  OOK  signal 
(Figure  17a)  show  that  the  modulated  carrier  ha.s  an  ON-OFF  characteristic.  The 
bit  rate  and  carrier  frequency  at  spectral  positions  5  and  11  can  be  estimated  using 
Figure  17b.  The  original  message  (within  a  180°  phase  uncertainty)  can  easily  be 
recovered. 

At  a  level  of  10-dB  (see  Figures  18a  and  18b)  the  noise  does  not  have  a  sig¬ 
nificant  effect  on  the  output  of  the  FFT  detector.  The  modulating  message  can  still 
be  recovered  (with  a  phase  uncertainty  as  mentioned  before).  But,  when  the  noise  is 
increased  to  a  level  of  0-dB  SNR  (Figure  19a  and  19b),  it  becomes  very  difficult  to 
distinguish  the  signal  from  the  noise. 

Figure  20  shows  the  output  averaged  over  frequency  and  time.  The  pres¬ 
ence  of  the  signal  can  be  recognized  at  all  noise  levels  and  the  carrier  frequency  and 
bandwidth  can  be  estimated  (see  la,  2a,  and  3a).  The  frequency  averaged  out¬ 
put  for  noise-free  OOK  (see  lb)  allows  the  message  code  to  be  estimated.  The 
amplitude  modulation  can  be  recognized,  and  the  bit  rate  can  be  obtained.  But, 
when  the  SNR  drops  to  10-dB  (see  2b),  it  is  difficult  to  estimate  OOK  features  from 
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Figure  16:  Averaged  output  of  FFT  detector  for  BPSK  signal  at  ail  SNR  levels. 
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Figure  20:  Averaging  output  of  FFT  detector  for  OOK  signal  at  all  SNR  levels. 
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the  averaged  output.  At  an  SNR  level  of  0-dB,  the  signal  vanishes  completely  (see 

3b). 

3.  FSK 

The  output  of  the  detector  driven  by  a  noise- free  FSK  signal  is  shown  in  Figure 
21a.  It  shows  a  signal  having  a  constant  amplitude  and  some  modulation.  Figure 
21b  shows  the  corresponding  contour  plot.  It  displays  FSK  characteristics  and  the 
bit  rate  can  be  extracted  from  this  plot  by  measuring  the  minimum  width  of  a  single 
message  bit.  The  carrier  frequency  and  the  frequency  shift  can  also  be  estimated  from 
the  contour  plot.  The  message,  within  a  sign  (±)  uncertainty,  can  be  recovered. 

Figure  22a  shows  the  output  for  the  FSK  signal  with  an  SNR  level  of  10-dB. 
Note  that,  even  at  this  level  of  noise,  we  can  still  recognize  the  presence  of  a  carrier 
and  we  note  that  the  amplitude  is  not  constant  due  to  the  added  noise.  Figure  22b 
allows  estimation  of  the  parameters  such  as  the  carrier  frequency,  the  frequency  shift, 
and  the  bit  rate.  It  still  permits  the  demodulation  of  the  message  (within  a  phase 
uncertainty  of  180“). 

Figures  23a  and  23b  shows  the  output  of  the  detector  driven  by  an  FSK  signal 
with  an  SNR  level  of  0-dB.  It  is  difficult  to  identify  the  FSK  signal.  The  output  can 
be  interpreted  as  OOK  at  a  carrier  frequency  of  5  and  13  (see  Figure  23b). 

Figure  24  shows  the  output  averaged  over  frequency  and  time.  The  time  av¬ 
eraged  output  (see  la,  2a,  and  3a)  shows  that  at  all  SNR  levels  a  non-symmetric 
spectrum  is  present,  which  is  probably  a  result  of  a  summation  of  more  than  one 
frequency  component.  One  of  the  carrier  frequencies  can  be  estimated  and  a  measure 
of  bandwidth  can  be  obtained.  Plot  lb,  2b,  and  3b  show  the  frequency  averaged 
output.  Plot  lb  (noise-free  case)  shows  a  constant  amphtude  carrier.  No  information 
can  be  obtained  at  10-dB  and  0-dB  (see  2b  and  3b). 
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Figure  23:  Output  of  FFT  detector  for  FSK  signal,  0-dB  SNR;  75%  overlap 
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Figure  24;  Averaged  output  of  FFT  detector  for  FSK  signal  at  all  SNR  levels 


4.  QPSK 

Figure  25a  shows  the  output  of  the  FFT  detector  driven  by  a  noise-free  QPSK 
signal.  It  shows  that  there  is  a  modulated  carrier.  Figure  25b  shows  the  corresponding 
contour  plot,  which  allows  bit  rate  estimation.  Comparing  this  figure  with  Figure  6, 
we  can  identify  a  phase  shift  of  180°  at  locations  19,  39,  47,  and  54,  and  phase  shifts 
of  plus  or  minus  90°  at  locations  2,  16,  24,  28,  32,  36,  44,  52,  and  58.  The  carrier 
frequency  at  positions  5  and  13  can  also  be  estimated. 

Figure  26a  shows  the  output  of  the  detector  at  an  SNR  level  of  10-dB.  The 
estimation  of  the  modulation  is  difficult  since  the  carrier  exhibits  the  characteristics 
of  FSK,  OOK,  or  BPSK.  Figure  26b  shows  the  corresponding  contour  plot.  The 
modulation  carrier  is  primarily  at  spectral  positions  5  and  13.  There  is  a  possible 
carrier  frequency  shift  at  some  locations  (i.e.,  at  delay  19  and  54).  Transition  points, 
probably  due  to  phase  shift,  are  present  and  can  be  used  to  estimate  the  bit  rate. 

Figures  27a  and  27b  show  the  output  of  the  detector  at  an  SNR  level  of  0-dB. 
It  is  impossible  to  identify  the  output  as  QPSK  or  any  other  modulated  signals  or  to 
estimate  any  parameters  due  to  the  noise. 

Figure  28  shows  the  output  averaged  over  frequency  and  time.  The  time  av¬ 
eraged  output  (see  la,  2a,  and  3a)  shows  that  the  signal  is  always  present  and  the 
carrier  frequency  and  bandwidth  can  be  estimated.  The  frequency  averaged  output 
in  the  noise- free  case  shows  constant  carrier  energy,  while  from  the  10-dB  and  0-dB 
SNR  plots  (see  2b  and  3b)  no  information  can  be  obtained. 

D.  TIME-FREQUENCY  DISTRIBUTION  METHODS 

1.  Method-1 

As  mentioned  earlier  (Eq.  4),  this  method  is  based  on  the  Fourier  transform 
of  lagged  products  (over  lag  m).  The  variable  n  is  an  independent  argument  here. 


42 


Figure  25:  Output  of  FFT  detector  for  QPSK  signal,  no  noise;  75%  overlap 


Figure  26:  Output  of  FFT  detector  for  QPSK  signal,  10-dB  SNR;  75%  overlap 
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Figure  28:  Averaged  output  of  FFT  detector  for  QPSK  signal  at  all  SNR  levels 


A  transform  length  of  16  is  assumed  and  x(n)  is  taken  as  the  sample  mean  over  these 
16  points.  The  time  axis  (samples  or  delay)  are  scaled  by  4,  corresponding  to  a  75% 
overlap  at  a  transform  length  of  16. 

The  outputs  of  the  method- 1  detector  are  shown  using  three  types  of  plots. 
A  three-dimensional  output  plot,  similar  to  a  spectrogram  of  the  data  shows  power 
versus  time  and  frequency  (Figures  29a,  30a,  31a,  33a,  34a,  35a,  37a,  38a,  39a,  41a, 
42a,  and  43a).  The  outputs  are  also  shown  as  contour  plots  for  more  analysis  and 
extraction  of  parameters  (Figures  29b,  30b,  31b,  33b,  34b,  35b,  37b,  38b,  39b,  41b, 
42b,  and  43b).  Also,  time  averaged  plots  (see  plots  la,  2a,  and  3a)  and  frequency 
averaged  plots  (see  plots  lb,  2b,  and  3b)  are  shown  (Figures  32,  36,  40,  and  44). 
a.  BPSK 

Figure  29a  shows  the  output  of  the  detector  for  a  BPSK  signal  with¬ 
out  noise.  Note  the  appearance  of  the  modulated  carrier  at  some  carrier  frequency. 
Fluctuation  in  the  magnitude  of  the  output,  may  indicate  the  presence  of  phase  mod¬ 
ulation.  The  contour  plot  in  Figure  29b  allows  the  estimation  of  the  carrier  frequency 
at  spectral  positions  5  and  13,  and  the  bit  rate. 

Figures  30a  and  30b  show  the  output  of  the  detector  for  a  signal  with  an 
SNR  of  10-dB.  Note  that  the  features  of  the  output  are  difficult  to  determine  in  both 
figures,  which  may  lead  to  an  improper  identification  of  the  type  of  modulation.  It 
is  obvious  that  the  output  of  the  detector  at  an  SNR  level  of  0-dB  (Figures  31a  and 
31b)  will  give  worse  results  than  those  shown  in  Figures  30a  and  30b. 

Figure  32  shows  the  output  averaged  over  frequency  and  time.  The  time 
averaged  output  (see  plot  la  for  noise-free,  2a  for  10-dB  SNR,  and  3a  for  0-dB  SNR) 
shows  that  at  each  of  the  different  noise  levels,  the  carrier  frequency  and  the  band¬ 
width  can  be  estimated.  The  frequency  averaged  output  shows  a  constant  output 
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Figure  29:  Output  of  Method-1  detector  for  BPSK  signal,  no  noise 


Figure  30:  Output  of  Method-1  detector  for  BPSK  signal.  10-dB  SNR. 


Figure  31:  Output  of  Method-1  detector  for  BPSK  signal,  0-dB  SNR. 
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for  the  noise-free  case  (see  lb).  The  output  at  an  SNR  of  10-dB  (2b)  and  an  SNR  of 
0-dB  (3b)  does  not  indicate  the  presence  of  a  signal. 

b.  OOK 

The  output  of  the  detector  for  the  OOK  signal  without  noise  is  shown  in 
Figures  33a  and  33b.  In  Figure  33a,  an  amplitude  modulated  carrier  is  present  and 
the  message  code  can  be  estimated.  Figure  33b  displays  the  corresponding  output  as 
a  contour  plot.  The  carrier  frequency  is  easily  verified  at  positions  5  and  13.  The  bit 
rate  can  be  estimated. 

Figure  34a  shows  the  detector  output  for  an  SNR  level  of  10-dB.  Note  that 
a  amplitude  modulation  can  be  observed  and  the  message  code  can  still  be  estimated. 
Figure  34b  allows  carrier  frequency  and  bit  rate  estimation. 

The  output  of  the  detector  at  0-dB  SNR  is  shown  in  Figures  35a  and  35b. 
Note  that  it  is  difficult  to  recognize  features  of  an  OOK  signal  due  to  the  high  noise 
level. 

Figure  36  shows  the  output  averaged  over  frequency  and  time.  The  time 
averaged  output  (see  plot  la  for  noise- free,  2a  for  10-dB  SNR,  and  3a  for  0-dB  SNR) 
shows  the  presence  of  the  signal,  and  the  carrier  frequency  as  well  as  the  bandwidth 
can  be  estimated.  The  frequency  averaged  output  is  shown  in  plots  lb  ,  2b,  and 
3b.  In  lb  the  output  is  noise-free  and  the  message  code  can  be  estimated  easily. 
The  amplitude  modulation  can  be  recognized  and  the  bit  rate  can  be  measured. 
In  2b,  where  the  output  has  an  SNR  of  10-dB,  OOK  type  modulation  can  still  be 
recognized  noting  instants  in  time  where  the  energy  is  minimal  (i.e.,  no  carrier).  The 
bit  rate  can  be  estimated  but  the  message  code  cannot  be  estimated.  The  frequency 
averaged  output  of  the  detector  (3b)  does  not  provide  modulation  clues  under  low 
SNR  conditions  (i.e.,  0-dB). 
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Figure  34:  Output  of  Method-1  detector  for  OOK  signal,  10-dB  SNR. 
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Figure  35:  Output  of  Method- 1  detector  for  OOK  signal,  0-dB  SNR, 


Amplitude  Amplitude 


SNR=lnfinity 


Frequency 

(la) 


SNR=10dB 


Frequency 

(2a) 


SNR=lnfinity 


0  20  40  60  80 

Samples 

(lb) 


SNR=10dB 


0  20  40  60  80 


Samples 

(2b) 


56 


c.  FSK 


The  output  for  the  noise-free  FSK  signal  is  shown  in  Figure  37a.  A  modu¬ 
lated  carrier  with  several  levels  of  magnitude  is  present.  Figure  37b  shows  the  carrier 
frequency,  shift  frequency  (positions  5  and  6,  12  and  13),  and  allows  bit  rate  estima¬ 
tion.  Figure  38a  shows  the  output  of  the  detector  for  an  SNR  of  10-dB.  A  modulated 
carrier  is  present.  Figure  38b  shows  the  FSK  feature  and  allows  carrier  frequency 
estimation  at  spectral  positions  5,  6,  12,  and  13.  The  bit  rate  can  still  be  estimated. 

Figures  39a  and  39b  show  the  output  for  an  SNR  level  of  0-dB.  We  cannot 
recognize  the  FSK  modulation  at  this  level  of  noise  and  this  may  lead  to  a  misclas- 
sification  of  the  modulation  (i.e.,  OOK  at  carrier  frequency  spectral  positions  5  and 
13). 

Figure  40  shows  the  output  averaged  over  frequency  and  time.  The  time 
averaged  output  (see  plot  la  for  noise-free,  2a  for  10-dB  SNR,  and  3a  for  0-dB  SNR) 
shows  that  at  the  different  SNR  considered,  the  carrier  can  be  recognized  having  a 
non-symmetrical  spectral  shape.  This  indicates  the  presence  of  more  than  one  spectral 
component,  which  is  a  feature  of  FSK.  A  measure  of  bandwidth  can  be  obtained.  The 
frequency  averaged  output  (see  lb,  2b,  and  3b)  shows,  in  the  noise-free  case  (lb),  the 
points  of  the  frequency  shift.  At  10  and  0-dB  SNR  even  that  clue  disappears  (see  2b 
for  10-dB  and  3b  for  0-dB). 
d.  QPSK 

Figure  41a  shows  the  output  of  the  detector  for  noise-free  QPSK  signal. 
It  appears  that  there  is  a  modulated  carrier.  Figure  41b  shows  that  the  bit  rate 
estimation  is  possible  and  that  transition  points  (i.e.,  phase  shifts)  can  be  identified 
as  a  180°  phase  shift  or  as  phase  shifts  of  90°.  The  carrier  frequency  can  be  estimated 
at  spectral  positions  5  and  13. 
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Figure  37:  Output  of  Method- 1  detector  for  FSK  signal,  no  noise 


Figure  38:  Output  of  Method- 1  detector  for  FSK  signal,  10-dB  SNR. 
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Figure  40:  Averaged  output  of  Method-1  detector  for  FSK  signal  at  all  SNR  levels 


Figure  41:  Output  of  Method- 1  detector  for  QPSK  signal,  no  noise. 


The  output  of  the  detector  for  an  SNR  of  10-dB  is  shown  in  Figures  42a 
and  42b.  A  modulated  carrier  centered  at  spectral  positions  5  and  13  is  present. 
Figure  42b  also  allows  an  OOK  interpretation  at  carrier  frequency  5  and  13. 

Figures  43a  and  43b  show  the  output  of  the  detector  for  an  SNR  of  0-dB, 
indicating  that  it  is  difficult  to  recognize  QPSK  at  this  low  SNR  level. 

Figure  44  shows  the  output  averaged  over  frequency  and  time.  The  time 
averaged  output  (see  plot  la  for  noise-free,  2a  for  10-dB,  and  3a  for  0-dB)  shows 
that  the  carrier  frequency  can  be  identified  and  the  bandwidth  can  be  estimated. 
The  frequency  averaged  output  is  shown  in  plots  lb,  2b,  and  3b.  For  the  noise-free 
signal  (see  lb)  a  constant  energy  is  shown.  For  10  and  0-dB  SNR  output  (see  2b,  3b, 
respectively),  the  signal  completely  vanishes  and  it  appears  that  only  noise  is  present. 

From  the  analysis  of  the  magnitude  of  the  method- 1  outputs,  we  note  that 
there  is  no  apparent  difference  in  performance  between  the  FFT  based  detector  and 
method- 1  even  though  some  estimated  parameters  differ  slightly. 

2.  Method-2 

As  mentioned  earlier  (Eq.  6),  this  method  is  related  to  WVD  and  is  based  on 
the  Fourier  transform  of  lagged  products  (over  time  n).  The  most  distracting  property 
of  WVD  is  that  it  can  have  spectral  components  in  regions  where  one  would  expect 
none.  These  spurious  values,  which  are  due  to  the  so-called  cross  terms,  are  particu¬ 
larly  prevalent  in  multi-component  signals.  Noise,  of  course,  is  the  superposition  of 
many  signals,  hence  causing  recognition  problems.  We  are  not  using  any  method  to 
minimize  or  filter  out  this  effect  at  this  time.  The  Choi-Williams  method  is  one  ap^ 
proach  to  deal  with  this  problem  [10].  This  problem  will  manifest  itself  in  the  detector 
output.  To  suppress  the  amplitude  belonging  to  the  negative  frequency,  complex  en¬ 
velopes  of  the  signal  are  used  in  the  program.  Note  that  the  reason  the  signal  appears 
at  twice  the  carrier  frequency  is  that  we  use  the  product  of  two  signals  in  the  time 
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Figure  42:  Output  of  Method-1  detector  for  QPSK  signal,  10-dB  SNR 


Figure  43:  Output  of  Method-1  detector  for  QPSK  signal,  0-dB  SNR 
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Figure  44:  Averaged  output  of  Method-1  detector  for  QPSK  signal  at  all  SNR 


domain  and,  as  expected,  in  the  frequency  domain  we  have  components  at  the  DC 
value  and  at  double  the  carrier  frequency  (i.e.,  sin^(a;f)  =  1/2[1  —  cos(2Ljt)]). 

The  output  of  this  detector  is  presented  using  two  types  of  plots.  First,  contour 
plots  are  shown  in  Figures  45  through  47,  49  through  51,  53  through  55,  and  57 
through  59.  Second,  outputs  averaged  over  time  (see  plots  la,  2a,  and  3a)  and  over 
frequency  (see  plots  lb,  2b,  and  3b),  are  given  in  Figures  48,  52,  56,  and  60.  Due 
to  the  selection  of  the  values  of  the  parameter  m,  the  frequency  averaged  outputs 
display  a  non-symmetry  and  the  number  of  computations  reduce  by  a  factor  of  2, 
with  delay,  in  the  range  —127  to  -M27  with  a  step  size  of  2. 
a.  BPSK 

Figure  45  shows  the  output  of  the  detector  for  a  BPSK  signal  without 
noise.  Note  the  presence  of  broadband  energy  between  spectral  positions  120  and 
140,  with  a  center  frequency  about  130.  A  symmetric  pattern  is  present  along  the 
delay  axis  between  the  two  limits  on  the  frequency  axis.  This  pattern  consists  of 
isolated  high  intensity  segments.  The  size  of  the  segment  is  not  equal  or  uniform, 
possibly  because  of  the  selection  size  of  variable  m.  Thus,  it  is  difficult  to  relate  this 
to  the  bit  rate  or  phase  shift. 

Figure  46  shows  the  output  of  the  detector  for  an  SNR  level  of  10-dB.  We 
note  that  the  presence  of  broadband  energy  between  the  approximate  limits.  The 
symmetric  pattern  is  still  present. 

Figure  47  shows  the  output  at  an  input  SNR  level  of  0-dB.  Note  that  some 
of  the  pattern  is  still  present  and  noise  and  cross  terms  are  now  dominant  about  the 
center  frequency. 

Figure  48  (see  la,  2a,  and  3a)  represents  the  time  averaged  output  for 
different  levels  of  SNR.  In  la  we  note  the  presence  of  a  carrier  frequency  and  spec¬ 
tral  spikes.  The  bandwidth  can  be  estimated.  The  presence  of  discrete  spectral 
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Frequency 


Figure  47:  Output  of  Method-2  detector  for  BPSK  signal,  0-dB  SNR. 
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Figure  48:  Averaged  output  of  Method-2  detector  for  BPSK  signal  at  all  SNR 
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components  may  be  due  to  the  cross  terms  of  the  WVD.  In  2a  we  get  a  similar 
result  as  la,  but,  in  3a  the  level  of  noise  masks  the  signal  completely.  The  frequency 
averaged  output  (see  lb,  2b,  and  3b)  provides  no  information  at  any  SNR  level. 

b.  OOK 

Figure  49  shows  the  output  of  the  detector  using  an  OOK  input  signal  un¬ 
der  noise-free  conditions.  Note  the  presence  of  broadband  energy  between  positions 
120  and  140  along  the  frequency  axis,  with  the  center  frequency  at  spectral  posi¬ 
tion  130.  A  symmetric  pattern  of  high  density  segments  of  energy  can  be  identified 
concentrated  at  the  center  frequency.  Patterns  of  this  form  may  indicate  OOK  type 
modulation. 

Figure  50  shows  the  output  for  an  SNR  level  at  10-dB.  Even  at  this  level 
of  noise,  patterns  similar  to  the  ones  obtained  in  Figure  49  can  be  obtained. 

Figure  51  shows  the  output  for  an  SNR  level  of  0-dB.  Note  that  at  this 
level  of  noise  most  of  the  signal  is  masked. 

Figure  52  shows  the  output  averaged  over  time  and  frequency.  The  time 
averaged  output  (see  la)  under  noise- free  condition  shows  a  strong  spectral  component 
at  location  130.  The  spectrum  of  the  signal  seems  to  be  noisy,  which  may  be  due 
to  cross  terms.  The  bandwidth  can  be  estimated.  At  an  SNR  of  10-dB  (see  2a), 
the  spectral  component  can  be  recognized  with  a  peak  at  the  appropriate  center 
frequency.  The  strength  of  the  cross  terms  are  more  evident.  At  an  0-dB  SNR  level 
there  is  no  trace  of  a  signal.  The  frequency  averaged  output  (see  lb,  2b,  and  3b) 
provide  no  useful  clues. 

c.  FSK 

Figure  53  shows  the  output  of  the  detector  for  an  FSK  signal  without 
noise.  Broadband  energy  is  present  between  spectral  positions  120  and  170.  Equi¬ 
distant  spectral  components  exist  with  spacing  equal  to  six  bins  on  the  frequency 
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Figure  49:  Output  of  Method-2  detector  for  OOK  signal,  no  noise. 
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Figure  52:  Averaged  output  of  Method-2  detector  for  OOK  signal  at  all  SNR  levels. 
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axis.  This  corresponds  to  the  bit  rate  of  the  modulating  signal.  Simulation  shows  a 
strong  sensitivity  of  these  spectral  peaks  to  the  number  of  transitions  in  the  original 
message  (i.e.,  modulation  index).  Also,  the  cross  terms  seem  to  be  present.  Segment 
shapes  between  spectral  positions  130  and  170  show  a  cosine  squared  behavior  along 
the  delay  axis.  Figure  54  shows  the  output  of  the  detector  an  SNR  level  of  10-dB. 
Results  similar  to  those  found  in  Figure  53  are  evident.  Figure  55  shows  the  output 
of  the  detector  at  0-dB  SNR.  There  are  still  some  high  energy  segments  of  spectral 
components  recognizable  along  the  delay  axis,  which  may  be  used  for  modulation 
identification. 

Figure  56  shows  the  output  averaged  over  time  (see  la,  2a,  and  3a)  and 
frequency  (lb,  2b,  and  3b)  at  different  SNRs.  At  noise-free  conditions,  the  time  aver¬ 
aged  output  (see  la)  shows  the  presence  of  spectral  components  with  constant  spacing 
corresponding  to  the  bit  rate.  At  an  SNR  of  10-dB  (2a),  the  spectral  components  can 
still  be  recognized  and  the  bit  rate  can  be  estimated.  At  an  SNR  level  of  0-dB  some 
of  the  spectral  components  are  still  present  and  the  bit  rate  can  still  be  estimated. 
The  averaged  frequency  plots  (see  lb,  2b,  and  3b)  convey  no  information  about  the 
modulation. 

d.  QPSK 

Figure  57  shows  the  output  of  the  detector  for  a  QPSK  signal  without 
noise.  Broadband  energy  between  spectral  positions  100  and  170  is  present.  The 
dominant  spectral  components  seem  to  be  centered  at  location  125.  A  symmetric 
pattern  along  both  axes  and  some  high  energy  spots  at  the  center  frequency  can  be 
recognized.  Figure  58  shows  the  output  for  a  signal  with  10-dB  SNR.  The  signal 
related  energy  is  concentrated  between  positions  120  and  140  and  the  noise  is  con¬ 
centrated  between  positions  70  and  180.  Figure  59  shows  the  detector  output  at  an 
SNR  level  of  0-dB.  The  features  of  the  signal  have  disappeared  completely. 
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Figure  56:  Averaged  output  of  Method-2  detector  for  FSK  signal  at  all  SNR  levels 


Figure  60  shows  the  output  averaged  over  time  and  frequency  at  different 
levels  of  SNR.  The  time  averaged  output  is  shown  in  la,  2a,  and  3a.  In  plot  la  (no 
noise  condition)  the  center  frequency  and  bandwidth  can  be  estimated.  The  presence 
of  spectral  spikes  can  be  recognized  and  may  be  related  to  cross  terms.  At  an  SNR 
level  of  10-dB  (see  2a)  similar  information  as  in  la  can  be  obtained.  At  an  SNR  level 
of  0-dB  the  detector  output  consists  mainly  of  noise.  The  frequency  averaged  plots 
convey  no  useable  information  (see  lb,  2b,  and  3b). 

The  method-2  output  permits  estimation  of  only  the  center  frequency  and 
the  bandwidth  (see  plots  averaged  over  time).  Method-2,  however,  provides  output 
patterns  that  are  unique  and  distinct  for  each  signal  modulation  at  all  SNRs  tested. 
This  may  lead  to  fingerprint  recognition  for  modulation  identification  and,  if  success¬ 
fully  implemented,  could  be  automated  to  work  in  conjunction  with  a  neural  network 
trained  to  recognize  (i.e.,  identify)  typical  modulation  types.  The  SNR  dependency 
can  be  improved  by  taking  longer  transforms,  provided  the  data  duration  permits  it. 
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Figure  60:  Averaged  output  of  Method- 
levels. 
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detector  for  QPSK  signal  at  all  SNR 
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V.  SUMMARY  AND  CONCLUSION 


The  signal  collection  capabilities  of  four  types  of  receivers  (interceptors)  con¬ 
sidered  in  this  work  are  compared  under  three  noise  conditions:  no  noise,  0-dB  SNR, 
and  10-dB  SNR.  Additionally,  they  are  compared  with  respect  to  their  response  to 
signal  modulation  and  their  dependence  on  detector  parameters  (e.g.,  FFT  size).  The 
four  receiver  types  are  designated  as  broadband  energy,  FFT  based,  method-1,  and 
method-2  detector. 

The  broadband  energy  detector  can  be  used  as  a  pulse-width  detector  for  OOK 
at  SNR  levels  greater  or  equal  to  10-dB. 

There  is  no  apparent  detection  capability  difference  between  the  FFT  based 
detector  and  the  method- 1  detector,  even  though  some  realizations  differ  slightly. 

The  application  of  method-2  can  be  equated  with  the  operation  of  a  circuit 
with  a  prefilter,  delay,  and  multiplication  as  shown  in  Figure  61.  In  some  sense,  this 
arrangement  is  similar  to  the  detector  advocated  in  [3].  One  of  the  main  differences 
is  due  to  the  FFT  taken  over  the  product  of  the  time  series  with  a  delayed  version  of 
itself. 

Table  1  shows  the  confusion  matrix.  It  shows  that  at  free-noise  conditions, 
all  modulations  can  be  identified  by  the  three  detectors.  At  a  level  of  10-dB  SNR, 
method-2  seems  to  be  more  powerful  than  the  other  two  detectors  in  identif5dng  the 
modulation.  The  FFT  fails  to  identify  QPSK  modulation  at  this  level  and  method-1 
fails  to  identify  BPSK  and  QPSK  modulation.  At  0-dB  SNR  all  detectors  fail  to 
identify  the  modulations.  Table  2  shows  which  of  the  parameters  can  be  extracted 
from  the  detectors  as  a  function  of  SNR.  The  carrier  frequency  can  be  identified  at 
all  SNR  levels  considered.  In  the  FSK  case,  only  one  of  the  carrier  frequencies  can  be 
estimated  and  identified.  The  FFT  based  detector  and  method-1  allow  estimation  of 
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most  parameters,  provided  the  SNR  is  sufficiently  large.  Method-2  is  more  useful  in 
estimating  the  type  of  modulation. 

Recommended  future  research  in  this  area: 

1.  More  than  one  data  realization  should  be  tested  to  get  statistical  information. 

2.  The  performance  of  the  detectors  should  be  investigated  more  thoroughly  by 
using  SNR  levels  in  smaller  steps  between  0-dB  and  20-dB. 

3.  All  simulations  are  done  with  detectors  synchronized  with  the  signals.  An  area 
for  future  research  is  to  evaluate  performance  of  the  detectors  having  different 
carrier  synchronization,  frequency  offset,  and/or  phase  offset. 

4.  All  detection  simulations  are  carried  out  on  a  signal,  in  addition  to  Gaussian 
noise.  An  area  for  further  research  is  the  extraction  of  parameters  from  one 
specific  signal  in  the  presence  of  an  additional  interfering  signal. 

5.  The  parameter  extraction  performance  of  method-2  should  be  investigated  after 
filtering,  that  is,  removal  of  cross  terms. 

6.  The  averaging  procedure  corresponds  to  the  DC  term  of  an  FFT.  The  time 
averaged  outputs  should  be  computed  by  taking  the  Fourier  transform  rather 
than  just  the  DC  term.  This  approach  may  reveal  modulation  parameters  as  a 
function  of  time. 
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Figure  61:  Method-2  block  diagram. 
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Table  1:  Confusion  Matrix 


Key  to  table: 


B 

O 

F 

Q 


BPSK  modulation  1 

OOK  modulation  * 

FSK  modulation 
QPSK  modulation  0 


identified 

multiple  modulation 

identified 

only  noise 


Note  :  An  empty  block  indicates  a  modulation  that 
cannot  be  recognized. 
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Table  2:  Parameters  Extracted  from  the  Modulations  after  Identification 


Key  to  table: 


fc 

carrier  frequency 

1  =  identified 

bandwidth 

*  =  multiple  modulation 

Br 

= 

bit  rate 

identified 

fsh 

frequency  shift 

0  =  only  noise 

Ps 

phase  shift 

Mp 

= 

message  pattern 
estimation 

Pu> 

pulse  width 

Note  :  An  empty  block  indicates  that  modulation 
parameters  cannot  be  recognized. 
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APPENDIX 


BPSK  Signal  Generation  Function  (bpsk.m) 

y.  Up  date  9.2.94 
y.  LCDR  ANWAR  AL-JOWDER. 

y.  This  function  generates  an  bpsk  signal  from  a  uniformly-distributed 
y,  random  sequence{-l,l}. 

y. 

cig 

clear 

N=input(' ENTER  THE  NUMBER  OF  SAMPLES ...  =  ') ; 
fs  =  input ('ENTER  THE  SAMPLING  FREQUENCY ; 


x=rand(l ,N) ; 

for  i=l:N; 

if  x(i)>.5 
x(i)=l; 

else 

x(i)=-l; 

end 

end 

y=boxcar(N) ; 
al=(y*x) ; 

al=reshape(al,l,(((N) .~2))); 

s=al . *sin((2*pi*f s/N) .*(1:((N).~2))); 


save  bpsk  s 


93 


OOK  Signal  Generation  Function  (ook.m) 

’/.  Up  date  9.2.94 

•/.  LCDR  ANWAR  AL-JOWDER. 

'/,  This  function  generates  an  bpsk  signal  from  a  uniformly-distributed 
y,  random  sequence{-l ,  1} . 

y. 

clg 

clear 

N=input(' ENTER  THE  NUMBER  OF  SAMPLES  ...  =  0 ; 
fs  =  input (CENTER  THE  SAMPLING  FREQUENCY ..= 0 ; 


x=rand(l ,N) ; 

for  i=l:N; 

if  x(i)>.5 
x(i)=l; 

else 

x(i)=0; 

end 

end 

y=boxcar(N) ; 
al=y*x; 

al=reshape(al,l, ((N) ."2)); 
s=al .*sin((2*pi*fs/N) . *(1 : ((N) . ~2))) ; 


save  ook  s  al 
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FSK  Signal  Generation  Function  (fsk.m) 


y,  Up  date  9.2.94 
y.  LCDR  ANWAR  AL-JQWDER. 

y.  This  function  generates  an  bpsk  signal  from  a  uniformly-distributed 
y,  random  sequence-C-1 , 1} . 

y. 


clear 

N=input(' ENTER  THE  NUMBER  OF  SAMPLES ...=') ; 
fs  =  input ('ENTER  THE  SAMPLING  FREQUENCY ; 
x=randn(l ,N) ; 
a=ones(l :N) ; 
for  i=l:N; 

if  x(i)>.5 
x(i)=l; 

else 

x(i)=0; 

end 

end 

y*boxcar(N) ; 

1=N~2; 

al=y*x; 

a2=reshape (al ,1,1); 
f=a2+f s; 

y,for  i=l  :length(al) ; 

s=sin(2*pi*(f ./length(a)) .*(1:1)) ; 

y.  end 

save  fsk  s 
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QPSK  Signal  Generation  Function  (QPSK.m) 

y.  Up  date  9.2.94 
y.  LCDR  ANWAR  AL-JOWDER. 

y.  This  function  generates  a  QPSK  signal  waveform.  Unif ormaly-distributed 
y.  random  sequence{-l ,  1} .  The  difference  between  the  four  levels  set  to 
y,  values  equal  to  {.25,  .5,  .75,  1.0}. 
clear 

N= input ( ' ENTER  THE  NUMBER  OF  SAMPLES ...='); 
fs  =  input (’ENTER  THE  SAMPLING  FREQUENCY ; 

y. 

b=N; 

x=rand(l ,b) ; 
for  i=l:b; 
if  x(i)<.25 

x(i)=2*pi; 

y. 

end 

end 

for  i=l:b; 
if  x(i)<.5 

x(i)=pi/2; 

y. 

end 

end 

for  i=l:b; 
if  x(i)<.75 

x(i)=pi; 

y. 

end 

end 

for  i=l:b; 
if  x(i)<=l,0 

x(i)=3*pi/2; 

y. 

end 

end 

y=boxccLr(b)  ; 
a2=y*x ; 

al=reshape(a2,l, ((b) .“2))  ; 
sl=sin((2*pi*f s/N) .*(0:length(al)-l)) ; 

s=sin(((2*pi*fs/N) .♦(0:length(al)-l))+al) ; 
save  qpsk  s  al 
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Broadband  Energy  Detector  (Energy.m) 

y.  Up  date  10.2.94 
y.  LCDR  ANWAR  AL-JOWDER 

y. 

y,  This  program  calculates  the  signal  energy , within  a  sliding  window, starting 
y,  from  any  sample  point  I  within  finite  length  of  256  samples. 

y. 

y. 

y. 

y,  List  of  variables/parameters  for  the  energy  program 

y,  U=  Numbers  of  the  window  samples  points  to  be  selected  by  the  user, 

y,  P=  Size  of  the  window  (in  sample  periods)  . 

y,  1=  Starting  Window  (in  sample  periods  relative  to  beginning  of  record)  . 

y,  YEN  =  Null  matrix. 

y. 

y. 

load  ook 
load  noiselO 
sn=s+W; 

y,  Set  the  input  maximum  SNR  in  dB. 

YEN=n; 

U=input(' ENTER  THE  SIZE  OF  N. . [8,16 ,32,64. .] ,U=  '); 

F=zeros(256,U) ; 
for  P=1:U; 

for  I=l:256-P+1; 

F(I,P)=0; 
for  i=I;I+P-l; 

F(I,P)=  F(I,P)+sn(i)'2; 

end 

end 

end 

YEN=[YEN;F] ; 

end 

save  energy  YEN 


FFT  Method  Detector  (shift2.m) 


*/.  Up  date  10.12.94 
*/.  LCDR  ANWAR  AL-JOWDER 

*/,  This  program  will  input  any  signal  generation  function  and  will  display 
*/,  the  FFT  detector  output  for  purposes  of  evaluations. 

•/. 

•/. 

*/,  List  of  variables/parameters  for  the  shift2  program 

'/,  G=  Null  matrix 

*/,  V=  Size  of  the  window  shift 

•/. 

clear 

y.cig 

load  bpsk 

’/.load  noise 
*/,sn=s+W ; 

*/,  The  parameter  V  is  the  size  of  the  FFT  window  selected  by  the  user 
y,  and  it  should  be  a  power  of  2  with  a  maximum  length  of  256. 

V=input( ' select  the  shift  point[l6,32, 128,256] ,V= 
y=[s  zeros  (1 :V)] ; 

G*C]; 

S=input(’Do  you  need  overlap. .?Y/N  CY]:’,'s'); 
if  S=='Y'; 

0L=input( 'Enter  the  Value  in  %  ,0L='); 
shift=((100-0L)/100)*(256/V) ; 

else 

shift=(256/V); 

end; 

for  i=l : shift : 256-V+l 

G=[G;  (fft(y(i:V+i-l)))  .''2]  ; 

end; 


save  shift  G 
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Correlation  Method  1  Detector  (Corrl.m) 

'/.  Up  date  10.2.94. 

•/.  LCDR  ANWAR  ALJOWDR. 

•/. 

'/,  This  prograjn  input  any  signal  generation  function 
'/,  and  will  display  the  spectral  correlation 

’/,  detector  output,  based  on  delay  time  n,  for  purposes  evaluation. 

•/. 


'/,  List  the  veriables/paxameters  for  the  corrl  program. 

'/,  Ycorl  =  Null  matrix  used  for  storing  the  data  signal. 

*/,  sn  =  The  cornibted  signal  with  noise. 

*/,  s  =  The  signal  without  noise. 

'/,  S2  =  The  correlated  signal  using  method  two  (delay  -  m)  . 

'/,  av  =  Simple  main. 

clear 

load  fsk 
'/.load  noiselO 
'/,sn=s+W; 

Ycorl=C]  ; 

S2=zeros(16) ; 
for  i*l:4:225; 

av=(s(i)+s(i+15))/16; 

S2=(abs(av*fft(s(i:i+15)))) .“2; 

Ycorl= [Ycorl ;S2] ; 


end 

save  corrl  Ycorl 
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Correlation  Method  2  Detector  (corr2.m) 

7.  Up  date  10.2.94 
•/.  LCDR  ANWAR  AL-JOWDER 

y. 

7,  This  program  will  input  any  of  signal  generation  function 
7,  and  will  display  the  spectral  correlation  detector  output 
7  for  purposes  of  evaluation. 

7,  This  equation  related  to  Wigner-Ville  Distribution. 

7. 

7.  List  the  veriables/parameters  for  the  Corrl  progreun. 

7,  Ycor2  =  Null  matrix  used  for  storing  the  signal  data. 

7.  sig  =  Input  signal  to-noise  power  ratio  in  dB. 

7,  xpl  =  Received  signal  (with  dely)plus  noise. 

7.  xp2  =  Received  signal  (without  delay)  plus  noise. 

7,  xp3  =  The  correlation  function  between  the  signal  and  the  delyed 

7.  signal. 

clear 
load  qpsk 

load  noiseO 
sn=s+W ; 

sl=hilbert (sn) ; 


Ycor2=  □  ; 
nn=mcLx(size(sl))  ; 

S1=C]; 

for  m=l-(255/2):2:(255/2)-l; 
if  m<0 ; 

xpl=Csl(l , l-m:nn)  zeros(l,-m)] ; 

else 

xpl=Czeros(l,m)  sl(l,l:nn-m)] ; 

end 

xp2=sl; 
xp3=xp2 .*xpl ; 


Sl=fft(xp3) ; 
Ycor2=[Ycor2;Sl] ; 
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end 


save  corr2  Ycor2 
x=linspace(-64,63,127) ; 
y=linspace(0,255,256)  ; 


101 


102 


LIST  OF  REFERENCES 


1.  Notes  for  EC  4590  (Digital  Satellite  Communication),  and  EC  2410  (Fourier 
Transformation),  Naval  Postgraduate  School,  1992  (unpublished). 

2.  Johnson,  V.  K.,  “Detection  of  Spread  Spectrum  Signals  in  the  Presence  of  Noise 
and  Narrowband  Interface,”  Master’s  Thesis,  Naval  Postgraduate  School,  Mon¬ 
terey,  CA,  Sep.  1986. 

3.  Kuehls,  J.  F.,  and  Evaggelos,  G.,  “Presence  Detection  of  Binary-Phase-Shift- 
Keyed  and  Direct-Sequence  Spread  Signals  Using  a  Prefilter  Delay  and  Multiply 
Device,”  IEEE  J.  on  Selected  Areas  in  Communication,  Vol.  8,  No.  5,  June 
1990. 

4.  Ziemer,  R.  E.  and  Peterson,  R.  L.,  Introduction  to  Digital  Communications, 
Macmillan,  New  York,  1992. 

5.  Roden,  M.  S.,  Digital  Communication  System  Design,  Prentice  Hall,  New  York, 
1988. 

6.  Spooner,  S.  G.,  “An  Energy  Analysis  of  the  Pseudo  Wigner-Ville  Distribution 
in  Support  of  Machinery  Monitoring  and  Diagnostics,”  Master’s  Thesis,  Naval 
Postgraduate  School,  Monterey,  CA,  June  1992. 

7.  Torrieri,  D.  J.,  Principles  of  Secure  Communications  Systems,  2nd  ed.,  Artech 
House,  New  York,  1992. 

8.  Proakis,  J.  G.,  Digital  Communications,  McGraw-Hill,  New  York,  1983. 

9.  Couch,  L.  W.,  II,  Digital  and  Communication  Systems,  4th  ed.,  Macmillan,  New 
York,  1993. 

10.  Cohen,  L.,  “Time- Frequency  Distributions  —  A  Review,”  Proc.  of  the  IEEE, 
Vol.  77,  No.  7,  July  1989. 


103 


104 


INITIAL  DISTRIBUTION  LIST 


No.  Copies 

1.  Defense  Information  Center  2 

Cameron  Station 

Alexandria,  VA  22304-6145 

2.  Library  Code  52  2 

Naval  Postgraduate  School 

Monterey,  CA  93943-5101 

3.  Chairman,  Code  EC/Mw  1 

Department  of  Electrical  and  Computer  Engineering 

Naval  Postgraduate  School 
Monterey,  CA  93943-5121 

4.  Professor  Ralph  Hippenstiel,  Code  EC/Hi  2 

Department  of  Electrical  and  Computer  Engineering 

Naval  Postgraduate  School 
Monterey,  CA  93943-5121 

5.  Professor  Donald  Wadsworth,  Code  EC/Wd  2 

Department  of  Electrical  and  Computer  Engineering 

Naval  Postgraduate  School 
Monterey,  CA  93943-5121 

6.  Professor  H.  H.  Loomis,  Jr.,  Code  EC/Lm  1 

Department  of  Electrical  and  Computer  Engineering 

Naval  Postgraduate  School 
Monterey,  CA  93943-5121 

7.  LTCOL  Yusuf  Mallalah  2 

Commander,  Bahrain  Amiry  Navy 

BDF,  P.O.  Box  245 
Naval  Base 
State  of  Bahrain 

8.  LTCOL  Abdullatif  Al-Zayani  1 

BDF,  P.O.  Box  245 

HQ,  Organization  &;  Planning  Department 
State  of  Bahrain 


105 


No.  Copies 


9.  LCDR  Anwar  Al-Jowder  2 

BDF,  P.O.  Box  245 
Naval  Base,  Electrical  Flotila  Office 
State  of  Bahrain 


106 


